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Imaging of the Bragg reflected x-ray beam is proposed and validated as an in-situ method for characteriza¬ 
tion of performance of double-crystal monochromators under the heat load of intense synchrotron radiation. A 
sequence of images is collected at different angular positions on the reflectivity curve of the second crystal and 
analyzed. The method provides rapid evaluation of the wavefront of the exit beam, which relates to local mis- 
orientation of the crystal planes along the beam footprint on the thermally distorted first crystal. The measured 
misorientation can be directly compared to results of finite element analysis. The imaging method offers an 
additional insight on the local intrinsic crystal quality over the footprint of the incident x-ray beam. 


I. INTRODUCTION 

Double-crystal x-ray monochromators utilizing nearly per¬ 
fect crystals of various semiconductor materials (primarily Si 
and Ge) are the primary optical devices for monochromati- 
zation of intense synchrotron radiation in the photon energy 
range of 5-30 keV. Such monochromators deliver monochro- 
matized radiation to a wide variety of x-ray experiments 
where relative photon energy bandwidths AE/E < 10“^ are 
required. Besides the energy bandwidth the main performance 
criteria of the double-crystal monochromators are the stabil¬ 
ity and reproducibility of the photon energy scale and beam 
position, monochromator throughput and spectral efficiency. 
Finally, an important criterion to address is the preservation 
of wavefront, which in the recent years has become increas¬ 
ingly important since coherence properties of x-ray sources 
are drastically improving. In this regard, in-situ characteri¬ 
zation and understanding of performance parameters of the 
double-crystal monochromators has become a high priority 
for the synchrotron community. 

A number of reports has been published over the years 
on the in-situ characterization of performance of various 
double-crystal monochromators. A major effort was related 
to minimization of thermal distortion on the first crystal 
of double-crystal monochromators for third-generation syn¬ 
chrotron sources (e.g., IH-Ill]). Solutions involving cryo- 
cooling of the monochromator crystals have been found and 
their performance limits investigated (e.g., |3-0])- To gain 
detailed understanding on the performance of high-heat-load 
double-crystal monochromators advanced characterization 
methods have been developed which combine finite-element 
analysis of thermo-mechanical response of the monochroma¬ 
tor crystals under the heat load of synchrotron radiation with 
the dynamical theory of x-ray diffraction for distorted crystals 

(ii). 

Advanced in-situ characterization techniques have been ap- 
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plied focused on the evaluation of angular profile of the ra¬ 
diation beam delivered by the monochromators. The actual 
methods that have been used included measurements of the 
monochromator crystal angular reflectivity curves (i.e., rock- 
ing curves) for different portions of the exit beam (e.g., lEB- 
O), measurements using crystal analyzers d, and x-ray 
grating interferometry . 

As commonly inferred, these techniques enable evaluation 
of the slope error on the first crystal due to the thermally in¬ 
duced distortion, which is a crucial step in the performance 
diagnostics of the high-heat-load monochromators. However, 
the demonstrated methods have not yet become widely spread 
in the community due to either a time consuming setup and 
measurement procedure or limited availability of specialized 
instrumentation such as crystal analyzers or x-ray interferom¬ 
eters. 

A commonly used basic method for in-situ performance 
characterization of double-crystal monochromators involves 
collection and interpretation of integrated (total) rocking 
curves at variable levels of power delivered by the x-ray 
sources. A typical experimental setup is schematically shown 
in Fig. |I1 The intensity of the exit beam dominated by the pri¬ 
mary low-order reflection (Si 111) is monitored using an ion¬ 
ization chamber ICO, while radiation reflected by the higher 
order reflection (e.g.. Si 333) is Altered and monitored using 
an additional ionization chamber ICl. The integrated intensity 
detected by the ionization chamber is plotted as a function of 
the Bragg angle of the second crystal which yields Si 111 and 
Si 333 rocking curves at various levels of power delivered by 
the x-ray source or power absorbed in the first crystal of the 
monochromator. 

In this work, an area detector (CCD) is placed after ICl to 
produce a sequence of images of the cross section of the exit 
beam (Si 333) at various angular positions of the second crys¬ 
tal on its rocking curve. It is shown how imaging of the exit 
beam yields additional insights on the in-situ performance of 
the monochromator. In particular, it is shown that (i) ther¬ 
mally stabilized crystals in the double-crystal monochroma¬ 
tors operate in the regime of weak lattice distortions where 
ray tracing approach is applicable, which permits unambigu¬ 
ous assignment of the portions of the images to the respec¬ 
tive portions of the crystal surface; (ii) at cryogenic temper- 
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FIG. 1: Experimental setup for in-situ performance characterization of double-crystal monochromators (see text for details). 


atures the angular profile of the exit beam in the scattering 
plane directly relates to the heat-load-induced misorientation 
of the crystal planes along the x-ray beam footprint; and (Hi) 
sequential imaging of the cross section of the exit beam yields 
topographs representing local misorientation of the lattice 
planes and local crystal quality over the incident beam foot¬ 
print. These conclusions are supported by experiments with 
double-crystal monochromators operated at a bending magnet 
beamline and an undulator beamline of the Advanced Pho¬ 
ton Source. The spatial resolution of the imaging method ap¬ 
proaches the x-ray extinction depth in Bragg diffraction from 
perfect crystals (typically a few microns). Empowered by 
diffraction imaging, rapid in-situ evaluation of the wavefront 
achieved in a single angular scan can facilitate design and im¬ 
plementation of adaptive monochromator systems where the 
state of the crystal is manipulated to compensate the induced 
thermal distortions. This can be considered as a viable ap¬ 
proach toward wavefront-preserving high-heat-load diffract¬ 
ing optics. 


II. THEORY 


The principle of operation of the monochromator relies on 
the double-crystal non-dispersive geometry shown in Fig.|2|a) 
where two identical Bragg reflections with equal d-spacing 
are utilized. If the two crystals are perfectly aligned (69 = 
0 ) the deflection of the incident x-ray beam with a finite di¬ 
vergence Ai/ji due to Bragg reflection from the first crystal is 
completely compensated by the deflection due to Bragg reflec¬ 
tion from the second crystal. The original direction of prop¬ 
agation centered about incident wavevector ko is preserved 
and the divergence of the incident x-ray beam is equal to the 
divergence of the exit beam (Aijji = A' 02 ). If an angular de¬ 
tuning (50 7 ^ 0 of the second crystal is introduced as shown in 
Fig. |3a) a partial overlap of the reflection regions of the two 
crystals takes place as shown in the cotTesponding DuMond 
diagram in the wavelength-angular exit space of the first crys¬ 
tal which is also the entrance wavelength-angular space of the 
second crystal (panel (rH-2)). The notations on the angular 
axes of the DuMond diagram is such that 9i — 9\ represents 
the entrance angular coordinate of the first crystal with respect 


to the center of its reflection region 0J. Similarly, O'l — 9\ is 
the exit angular coordinate of the first crystal, 02 — 6*2 the 
entrance angular coordinate of the second crystal and 9'^ — 9^ 
is the exit angular coordinate of the second crystal. Since the 
second crystal is rotated clockwise 02 = 0o + (50. The wave¬ 
length axis of the DuMond diagrams is given in the relative 
units AA/Ao = (A — Ao)/Ao. 

As a result of the partial overlap, the center of the reflection 
region in the exit space of the second crystal (panel ( 2 ’)) shifts 
by (50 and the central portion of the exit beam (denoted by 
wavevector k^) propagates at an angular offset 250 from ko 
as shown in the ray tracing diagram of Fig.|2ja). 

The width of the rocking curve of the second crystal results 
from the convolution of the two identical reflection regions: 


A 02 = '/2A9}{, A9h = e/f tan0o (1) 

In this expression A9h is the intrinsic angular width of the 
Bragg reflection (also known as Darwin width), eh is the in¬ 
trinsic relative energy width of the Bragg reflection (a quantity 
given by the dynamical theory of x-ray diffraction), and 0 o is 
the central glancing angle of incidence related to the center Aq 
of the selected wavelength region via the Bragg’s law. 


Ao = 2dH sin 00 , 


( 2 ) 


where dn is the d-spacing of the working reflection. 

In the practical case when the relative photon energy band¬ 
width of the x-ray source is greater than eh, the relative pho- 
ton energy bandwidth selected by the monochromator is IBlIl 


|AA| _ Atpi 
Ao tan 00 


(3) 


This relationship can be easily illustrated using the Du¬ 
Mond diagram in panel (1) of Fig. Wia). The total relative 
bandwidth is the combination of the intrinsic relative band¬ 
width Eh and the divergence of the incident x-rays Aijji ex¬ 
pressed in relative energy units using the linearized Bragg’s 
law. 

X-ray absorption of the intense incident x-ray beam in the 
first crystal and dissipation of the absorbed energy leads to a 
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FIG. 2: Ray tracing and DuMond diagrams of the double-crystal monochromator (see text for more details): 

(a) The case of two perfect crystals with an identical symmetric crystallographic orientation. If an angular offset 59 is introduced a partial 
overlap of the reflection regions of the two crystals takes place in (I’-l-T). 

In cases (b) and (c) the distorted surface of the first crystal is considered as a superposition of flat perfect segments with variable angular 
orientation. X-rays reflected from these segments are schematically shown by variable color. 

(b) The case of low order reflection (Si 111). The individual segments are presented as a sequence of strongly overlapping reflection regions. 

(c) The case of high order reflection (Si 333). For simplicity the entrance space of the first crystal is presented as a superposition of non¬ 
overlapping local reflection regions. The reflection region of the second crystal selects a portion of radiation which corresponds to a segment 
with a particular local misorientation depending on the angular offset 59 . 
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non-uniform temperature distribution and thermal distortion 
of the crystal (known as the ’’heat bump”) shown schemati¬ 
cally in Fig.|2jb,c). In the following we will show that under 
certain conditions the image of the exit beam cross section 
is directly related to the thermal distortion of the first crystal 
along the footprint of the incident beam. 

The non-uniform temperature distribution and the result¬ 
ing thermal distortion lead to a variation in the d-spacing and 
in the local orientation of the crystal planes involved in the 
diffraction process. Differentiation of the Bragg’s law yields; 

5A Sdn 59 
Aq (Ih tan 00 

At each fixed photon energy a variation in the d-spacing 
Sdn yields an angular deviation 

(50 =tan 00 (5) 

dn 

The influence of such angular deviation on the resulting 
wavefront distortion can be neglected at cryogenic tempera¬ 
tures in vicinity of zero thermal expansion for the crystal ma¬ 
terial. This approximation is quite reasonable based on the 
fact that thermal expansion is rather small at cryogenic tem¬ 
peratures and an assumption that temperatures of the working 
crystal region do not deviate substantially from the zero ther¬ 
mal expansion temperature. To illustrate this approximation, 
we consider a Si crystal of a high-heat-load monochromator 
under typical operating conditions at a third-generation syn¬ 
chrotron source. We assume that the maximum temperature 
variation on the hrst crystal along the beam footprint is on 
the order of 10 K in vicinity of the Si zero thermal expan¬ 
sion point (Tq Ri 124 K), and that the operating photon energy 
Eq = hc/Xo — 5 keV (tan0o ~ 0.4). This brings the crystal 
close to the upper limit in the operational heat load (r; 400 W 
of absorbed power delivered by undulator A at the Advanced 
Photon Source, undulator deflection parameter K r 2). Tak¬ 
ing into account the temperature dependence of the Si lattice 
parameter ifl^ Eq. |5]yields (50 < 1 prad, which is small com¬ 
pared to the characteristic divergence of the synchrotron un¬ 
dulator source (r 10 /rrad). We note that in a general case the 
results of measurements in the double-crystal geometry are si¬ 
multaneously affected by the variations in the d-spacing and 
the misorientation of the reflecting planes. Thus, the mea¬ 
sured angular profile represents a gradient in the wavefront 
of the outgoing beam. Deconvolution of the two contribu¬ 
tions can be accomplished using monochromatization of the 
incident beam and application of methods of high-resolution 
x-ray diffraction topography (e.g., ifl^ l. This however, is not 
practical for the in-situ performance diagnostics. 

Neglecting the angular deviation due to temperature- 
induced variation in d-spacing for each photon energy implies 
that the wavefront is primarily affected by the local misori¬ 
entation of the Bragg planes. Here, we need to make another 
important assumption that the crystal is only weakly distorted, 
such that geometrical optics (i.e., ray-tracing) approach is ap¬ 
plicable. Strong crystal distortions lead to generation of new 
wavefields due to interbranch scattering (e.g, |IMII3,II11])- As 


TABLE I: Numerical values of the right-hand side and the left-hand 
side of Eq. [^evaluated separately for the Si 333 reflection at different 
photon energies to illustrate applicability of the weakly distorted 
crystal approximation. 

= hc/Xo - photon energy, 

Aq - extinction length, 

^ A 0333 - angular half-width. 

[keV] 15 21 27 33 39 

[prad] 0.1 0.07 0.06 0.05 0.04 
iA0333 [/rrad] 1.7 1.2 0.9 0.7 0.6 


a result it may become difficult to draw a cotTespondence be¬ 
tween a local region in the cross-section of the exit beam and 
a local region on the crystal surface. The physical interpre¬ 
tation of the criterion for the smallness of crystal distortion 
where the ray-tracing approach is applicable is given by Au- 
thier ifT^ . The lattice misorientation over the extinction length 
should be much smaller than the angular half-width of the in¬ 
trinsic Bragg reflectivity curve; 


9(A0) 

dl 


Ao <C 2 ^^^' 


(6) 


The extinction length is a minimal characteristic length from 
which x-rays are dynamically reflected. In the Bragg geome¬ 
try it is given by the real part of the following expression. 


Aoi/7o|7/»I 

\C\VXhX-h’ 


(7) 


where 70 and 7 /j are the direction cosines of the incident and 
reflected waves respectively, Xh and x-h are the Fourier com¬ 
ponents of dielectric susceptibility and C is the polarization 
factor (C=l for a-polarization of the incident radiation). The 
intrinsic angular width A9h as well as the extinction length 
Aq are characteristic to a chosen working reflection. 

As frequently reported in the literature (e.g., sain [II), 
for the cryogenically cooled Si monochromators the total mis¬ 
orientation of the Bragg planes along the beam footprint on 
the first crystal is 0 r 10 /rrad. The variation in the lattice 
misorientation in Eq. | 6 ]is on the order of 0 sin Oq/A y, where 
Ay R 1 mm is the typical size of the incident x-ray beam in 
the scattering plane. Table [B] shows values of the right-hand 
side of Eq.| 6 ]evaluated for the most practical Si 333 reflection 
at different photon energies. The values of the angular half¬ 
width of the reflection A A 0333 are also given for comparison. 
Thus, the criterion of weakly distorted crystal is satisfied for 
the high-heat-load monochromators under the operating con¬ 
ditions. 

Within the framework of the ray-tracing approach the ther¬ 
mally distorted first crystal can be considered as a superpo¬ 
sition of flat perfect segments with a variable angular orien¬ 
tation. X-rays reflected from these segments are schemati¬ 
cally shown by variable color in the ray-tracing diagrams (Fig. 
|2b,c)) and in the corresponding DuMond diagrams as a se¬ 
quence of overlapping reflection regions. The second crystal 
in the double crystal scheme is assumed to be perfect since 
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it is subjected to a substantially smaller heat load of the x- 
ray beam Bragg-reflected from the first crystal and scattered 
radiation. Thus, the reflection region of the second crystal re¬ 
mains unchanged (grey shaded area in DuMond diagrams of 
Fig. El). Rotation of the second crystal results in the shift of 
its reflection region across the overlapping regions represent¬ 
ing portions of radiation emanating from the first crystal (pan¬ 
els (l’-i-2) in Fig.|2lb,c)). The resulting interaction is drasti¬ 
cally different for the case of low order reflection (e.g.. Si 111) 
shown in Fig. Elb) and the case of higher order reflection (Si 
333) as shown in Fig. Etc). These two interactions occur si¬ 
multaneously at different photon energies and can be observed 
separately as discussed in the previous section. 

The variation in the misorientation of the first crystal in¬ 
teracting with the incident x-ray beam can be described by 
a maximum peak-to-valley deviation 0 from the nominal 
’’zero” value that corresponds to the perfect flat crystal ori¬ 
ented at a glancing angle of incidence 9q with respect to the 
incident beam. 

The intrinsic angular acceptance of the low-order reflection 
> 0. The reflection regions corresponding to different 
segments of the crystal strongly overlap in the wavelength- 
angular entrance space of the first crystal as shown in panel 
(1’) of Fig. |2jb) (only the boundaries of the individual re¬ 
gions are shown by varying color; the overlap region is shaded 
by cyan color). Without loss of generality it is assumed that 
the misorientation is a monotonic function of a spatial coor¬ 
dinate along the footprint of the incident beam. The diver¬ 
gence of x-rays reaching each segment of the entrance sur¬ 
face of the first crystal is small (within the framework of 
geometric optics it approaches zero with the size of the seg¬ 
ment). The wavelength-angular space selected by the total 
incident divergence Ai/;i can be approximated by a parallelo¬ 
gram ABCD where the points A, B, C, and D are defined as 
intersections of the angular boundaries limiting the total di¬ 
vergence with the boundaries of the reflection regions of the 
furthermost segments. In this simplified consideration the to¬ 
tal area of ABCD scales as ((0 + Ai/ji)/ tan Oq + eh)- Thus, 
a distorted crystal will select a greater photon bandwidth from 
the incident beam as compared to the bandwidth selected by 
the perfect crystal (Eq. El>- The illuminated total reflection 
region ABCD is sliced and consequently projected onto the 
wavelength-angular space (1 ’ -i-2) for each individual segment. 
For simplicity zero angular displacement of the second crystal 
from the central angle 9o is assumed {69 — 0). Upon sub¬ 
sequent projection onto the exit space of the second crystal 
(panel (2’)) it is seen that the reflected intensities correspond¬ 
ing to the different segments are present simultaneously in the 
reflected beam and propagate in different (divergent) direc¬ 
tions. The angular divergence of the exit x-ray beam is given 
by 

AiIj 2 = A'tpi + 0 ( 8 ) 

This increase in the angular divergence is critical for eval¬ 
uation of performance of synchrotron beamlines and is often 
interpreted as a virtual source located at a closer distance. In 
the context of the double-crystal geometry the misorientation 
of the Bragg planes on the first crystal due to thermal distor¬ 


tion represents slope error on the surface of the first crystal 
and is the origin of the additional divergence. The spatial dis¬ 
tribution of the slope error variation along the beam footprint 
can be extracted from the low-order reflectivity if the measure¬ 
ment method has angular selectivity. This is accomplished in 
prac tice using crystal analyzers ifl^ or grating interferometry 

In the case of a high-order reflection (e.g.. Si 333) A9^^^ < 
0 a qualitatively different picture emerges (Fig. Elc)). For 
simplicity it is assumed that the entrance space of the first 
crystal for the higher order reflection can be represented as 
a superposition of non-overlapping local reflection regions. 
These regions are sequentially illuminated with the divergent 
fan of the incident radiation. The line AB shown in panel (1) 
serves as a guide to the eye for subsequent illumination of the 
non-overlapping crystal regions. It can be seen that the region 
ABCD of Fig. Etb) transforms into the line AB under the as¬ 
sumption of zero overlap between the regions and the fraction 
of the incident divergence A'lpi (corresponding to an individ¬ 
ual crystal segment) approaching zero. In the exit space of the 
first crystal (l’-i-2) the entrance reflection region of the sec¬ 
ond crystal selects a portion of radiation which corresponds 
to a certain illuminated region representing a particular lo¬ 
cal misorientation (or a small range of misorientations) on the 
first crystal depending on the angular offset 69. The resulting 
reflection region projected onto the exit space of the second 
crystal (panel (2’)) has an angular width 

AV^2 = A6|333 (9) 

This selected portion of radiation propagates with an angu¬ 
lar offset 69 from the original direction ko as shown in the 
ray tracing diagram of Fig. EJc). We note that Eq. |9] is an 
idealization related to the assumption of non-overlapping dis¬ 
crete regions. In reality, the exit divergence will depend on 
the relationship of the local slope error to the intrinsic angular 
acceptance of the high-order reflection A9^^^. 

Eor the purpose of diagnostics using x-ray diffraction imag¬ 
ing the choice of a high order reflection such as Si 333 is par¬ 
ticularly convenient since its intrinsic angular width {A9^^^ ~ 
0.5 — 3 /rrad) is smaller than the expected total variation in 
misorientation 0 while the extinction depth is not particularly 
large, which ensures that the criterion given by Eq. |6]is satis¬ 
fied. 

Erom the analysis given above it is clear that in the absence 
of plane misorientation all illuminated regions of the first crys¬ 
tal will equally contribute to the image of the exit beam-cross 
section in the observation plane at any angular offset 69 of the 
second crystal. In contrast, if thermal distortion is present re¬ 
gions corresponding to different misoriented segments of the 
first crystal will be highlighted by maximum intensity at dif¬ 
ferent angular positions of the second crystal 69. In the case 
of a monotonic distribution of misorientation as a function of 
position on the first crystal along the beam footprint the maxi¬ 
mum intensity in the resulting image will monotonically shift 
in the vertical direction as 69 is scanned. 

To draw one-to-one correspondence between a particular 
location on the crystal and the location in the observation 
plane the observation plane should be located at a distance 
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D < Sx/A-tpi where Sx is the required spatial resolution. Al¬ 
ternatively, one can scale the image in the observation plane 
by the factor L/{L + D) to obtain the images of the crystal. 
It is also noted that the field of view is typically limited by an 
aperture of size Ay placed upstream the monochromator. In 
order to avoid diffraction on the aperture it is desirable to limit 
the observation distance D <C Aj/^/Aq. For typical aperture 
sizes Ay Ki 1 mm and Aq ~ lA aperture-related diffraction 
effects are expected for D Ri 10^ m. 

Quantitative information on the state of the first crystal can 
be extracted from the sequence of images obtained at different 
56 using rocking curve imaging strategy (e.g., lEoLI^ I. The 
images are sorted on per-pixel basis and local rocking curves 
are analyzed for each pixel. Rocking curve topographs can be 
constructed as maps of local reflected intensity, curve’s width 
and peak position. The map of the reflected intensity repre¬ 
sents the beam footprint on the first crystal superimposed on 
the local Bragg reflectivity. The map of peak positions rep¬ 
resents gradient of the wavefront in the scattering plane (i.e., 
the tangential direction) (H]. The map of the width provides 
an insight on crystal quality. For slightly distorted perfect 
crystals the expected values should be close to the theoreti¬ 
cal value In the following analysis of angular and 

spatial resolution in the rocking curve topographs is presented. 


III. ANGULAR SENSITIVITY AND SPATIAL 
RESOLUTION 

In the double-crystal configuration the second crystal plays 
the role of a Si 333 analyzer, which selects different portions 
of the x-ray fan emana ting from the first crystal as also noted 
in earlier studies (e.g., fl^ ). In principle, the sensitivity of an¬ 
gular measurements using the Si 333 analyzer is not directly 
related to the intrinsic angular width of the reflection A9^^^. 
It is rather defined by the measurement accuracy of peak po¬ 
sitions of the local rocking curves. This accuracy generally 
depends on the signal statistics and precision/stability of the 
crystal rotation. In the case of double-crystal monochromator 
the signal statistics is excellent due to the high photon flux of 
the exit beam. The angular stability may approach the level 
of 0.1 /irad which is about the practical limit for the angular 
sensitivity. 

The spatial resolution is limited by a number of geometrical 
factors. In practice the diffraction limit can be excluded from 
the consideration since it is small compared to another funda¬ 
mental quantity, the wave held penetration depth into the crys¬ 
tal. For an incident wave in Bragg diffraction the penetration 
depth approaches the extinction depth which attains the min¬ 
imal value Ao/27r in the center of the reflectivity curve. This 
sets a fundamental limit in the spatial resolution as shown in 
Fig. [3] On the tails of the reflectivity curve the penetration 
depth into the crystal approaches the absorption depth C, (typ¬ 
ically C > Ao/27r for practical cases in hard x-ray optics). 
Thus, in the measurements of peak position of local rocking 
curves the spatial resolution in the observation plane due to 


FIG. 3: Ray tracing diagram in Bragg diffraction illustrating spatial 
resolution due to finite penetration depth of the wave field into the 
crystal. 



the finite penetration depth approaches 

=-Aocos6»o (10) 

TT 

For Si, the primary crystal used in high-heat-load 
monochromators the extinction depth of the 333 reflection is 
about 4 ym. In the practical range of the primary photon en¬ 
ergies 3-15 keV, cos 00 — 1 and 5yp ~ 4 pm. 

Another source of blurring in the images of the exit beam 
is the finite source size S. The related spatial resolution in the 
observation plane at a distance D from the monochromator is 

5ys = D^. ( 11 ) 

For an undulator source at a third generation synchrotron 
S' ~ 25 pm (FWHM) while for a bending magnet beamline 
S ~ 70 pm (FWHM). 

Also, the intrinsic angular deviation of the exit beam, which 
is present locally for each segment and illustrated for perfect 
crystals in Fig. |2ja) has to be taken into account. The related 
spatial resolution in the observation plane is 

5ya = D2V2A0^^^. (12) 

The contributions to spatial resolution given by Eqs. (dUl 
and (fT^ can be optimized by placing the area detector at a 
sufficiently small distance from the monochromator D L. 
The resulting spatial resolution can be estimated as 

Sy= ^ 5yl + 5yl + 5yl. (13) 

Finally, it is noted that diffraction from single crystals leads 
to spatial Altering of the object wave ll2^ . Precise interpreta¬ 
tion of crystal imperfections in diffraction im aging ge nerally 
requires solutions of Takagi-Taupin equations UM . 

IV. EXPERIMENTAL 

Sequential diffraction imaging of the exit beam of the high- 
heat-load double-crystal Si 111 monochromators was per¬ 
formed at the Advanced Photon Source using experimental 
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TABLE II: Experimental parameters during x-ray diffraction imag¬ 
ing experiments of double-crystal high-heat-load monochromators at 
1-BM and 20-ID beamlines of the Advanced Photon Source: 

L - distance from the source to the monochromator, 

D - distance from the monochromator to the observation plane, 

- photon energy of the working Si 111 reflection, 

P - total radiation power absorbed in the Si crystal, 

Sy - spatial resolution in the tangential direction (Eq. in}, 

T - measured reference temperature on the first monochromator crys¬ 
tal. 

Beamline L [m] D [m] [keV] P [W] 6y [/rm] T [K] 
1-BM 27.5 7.0 6.0 1.3 58 289 

20-ID 30.0 20.0 11.0 93.5 84 82 

setup shown in in Fig. |I] The exit beam was filtered to se¬ 
lect Si 333 reflection for imaging. In the first experiment 
a water-cooled monochromator operating at a bending mag¬ 
net 1-BM beamline was studied. The working Si 111 re¬ 
flection of the monochromator was set to select photon en¬ 
ergies =6 keV from the incident broadband spectrum 
of bending magnet radiation. The cross section of the inci¬ 
dent beam was limited with 1x1 mm^ aperture located up¬ 
stream of the monochromator at a distance of 23 m from the 
source.Stabilization was ensured by repetitive measurements 
of the Si 333 rocking curve using an integrating detector until 
narrow curves were obtained with FWHM being close-to the 
theoretical value of 4.2 /rrad. 

In the second experiment a cryo-cooled (indirect liquid ni¬ 
trogen cooling) double-crystal monochromatorat 20-ID undu- 
lator beamline was studied. The working Si 111 reflection of 
the monochromator was set to a photon energy Eq^^ =11 keV 
to select the primary harmonic of the undulator that was 
tuned accordingly (undulator type A at the Advanced Photon 
Source). 

In both experiments prior to imaging sufficient time was 
allowed for thermal stabilization under the heat load of the in¬ 
cident x-rays. A sequence of images of the Si 333 reflection 
was collected corresponding to different angular positions of 
the second crystal on its rocking curve. The temperatures of 
the crystals were monitored using thermocouples directly af¬ 
fixed to the crystal base in each monochromator. The power 
absorbed by the first Si crystal P was calculated using mass- 
energy x-ray absorption coefficient of Si and equations for 
synchrotron radiation 125] with the relevant parameters 
of a bending magnet and undulator A at the Advanced Pho¬ 
ton Source. The main experimental parameters are listed in 
Table llVI for each experiment. 

V. ANALYSIS AND DISCUSSION 

The collected images were sorted to calculate local rock¬ 
ing curves for each pixel. Sequential x-ray topo grap hs were 
computed using rctopo code of DTXRD package 112'^ . To re¬ 
duce noise in the topographs the parameters of rocking curves 
at each pixel were extracted using fitting with the Gaussian 
function. Figures |4|a-c) show sequential x-ray topographs 
of the first experiment with the water-cooled double-crystal 


monochromator at the bending magnet beamline. Each to¬ 
pograph is sliced as shown by the vertical dashed line. The 
resulting distributions are plotted in Figs. Hfd-f) accordingly. 
Figures |4|a,d) show distributions of the normalized peak in¬ 
tensity of local rocking curves. Aside from a few local dis¬ 
turbances the peak intensity is quite uniform as expected for 
high-quality Si monochromator crystals at low heat load. Dis¬ 
tributions of the rocking curve width shown in Fig. |4|b,e) 
reveals more detailed information on the local disturbances 
suggesting presence of intrinsic localized strain in several lo¬ 
cations (increased local FWHM as compared to the theoret¬ 
ical value of A02 — 4 prad for Si 333 double-crystal rock¬ 
ing curve at 18 keV). For small regions of interest such as the 
one defined by the entrance aperture these disturbances can be 
avoided using crystal translations. Distributions of the rock¬ 
ing curve peak position shown in Fig. life,!) reveal a small 
gradient of ~ lp,rad peak-to-valley. We note that in the case 
of water-cooled monochromator crystal the thermal expansion 
coefficient of Si near the room temperature is quite substantial 
a ~ 2.6 X 10“®. Small temperature variations along the beam 
footprint (e.g., ST « 1 K) can lead to a substantial variation 
in the d-spacing Sdn/dn = ctST r; 10“®. Therefore, the an¬ 
gular deviation defined by Eq. |5]can no longer be neglected. 
Thus, at room temperatures the measured gradient can not be 
solely attributed to the misorientation of the Bragg planes. 
Nevertheless, regardless of the origin the measured gradient 
represents the disturbance in the wavefront due to non-ideal 
state of the crystal. 

Eigures |3a-c) show sequential x-ray topographs of the 
second experiment with the cryo-cooled double-crystal 
monochromator at the undulator beamline. Similarly to the 
previous case each topograph is sliced as shown by the ver¬ 
tical dashed line. The resulting distributions are plotted in 
Eigs. |5ld-f) accordingly. The topographs show a portion of 
the imaging plane which includes the central radiation cone. 
The map of the rocking curve peak intensity Eigures |3 a) rep¬ 
resents the envelope of the Si 333 reflected intensity. At each 
particular angle on the rocking curve the reflected intensity is 
present only in a narrow stripe in the imaging plane as shown 
by the sequence of alternating images (Supplementary sec¬ 
tion). This stripe moves in the vertical direction as angle of 
the second crystal is scanned. This behavior manifests itself 
as a strong gradient shown in Eigs. HJc,!) (« 8 prad peak- 
to-valley). The width of the local rocking curves is quite 
uniform (Eigures Sb, e)) and close to the theoretical value 
A02 = 2.2 prad. This indicates high crystal quality over the 
entire observation region. 

As discussed in section |II] the measured angular profile of 
the rocking curve peak position should represent the distribu¬ 
tion of heat-load-induced misorientation of the Bragg planes 
or the surface slope error given that the initial state of the crys¬ 
tal in the absence of heat load is strain-free. It is imperative 
to compare the measured slope error with predictions of finite 
element analysis (EEA). A volumetric heat source in the Si 
crystal produced by the undulator beam was obtained using 
equations describing spectral/spatial distribution of undulator 
radiation ll^ and the energy-dependent x-ray absorption co¬ 
efficient of Si. The heat source was introduced in finite el- 
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FIG. 4: Rocking curve topographs of the exit beam of a double-crystal monochromator (tuned to a photon energy of = 6 keV) at a 
bending magnet beamline of the Advanced Photon Source (1-BM): a map of the normalized peak intensity of the local rocking curves (a), 
a map of the curve width (b) and a map of the peak position (c). Each map is sliced as shown by the vertical dashed line. The resulting 
corresponding distributions along these lines are those of peak intensity (d), curve width (e) and peak position (f). The theoretical value of 
the rocking curve width for perfect crystals in the double-crystal geometry (4 /irad) is shown by the dashed line in (e). The dashed line in (f) 
represents the reference zero level. The size of the incident x-ray beam is limited by 1x1 mm^ aperture upstream of the monochromator. 


ement analysis with boundary conditions consistent with the 
typical scheme of an indirect-cooled Si monochromator ||28|] . 
The cooling channel of the monochromator crystal holder had 
wire-coil inserts for enhanced heat transfer such that the ef¬ 
fective transfer coefficient was assumed to be 8 kW/m^K“^ 
. The heat transfer in the silicon-copper interface with in¬ 
dium foil as the interstitial material is known to depend on the 
pressure applied ifsol] . The effective value of 4 W/m^K“^ was 
assumed for this parameter in the performed finite element 
analysis. The slope error along the center line of the beam 
footprint on the first crystal extracted from the results of FEA 
is shown in Fig.|6]by the green dashed lines. 

The measured angular profile of the central slice (Fig¬ 
ure Elf)) was recalculated in the coordinate along the sur¬ 
face of the crystal (1). The profile was aligned with respect 
to the reference zero level of the central radiation cone using 
the bottom edge of the aperture, which was clipping the x-ray 
beam incident on the monochromator. The angular profile is 
shown in Fig. |6ta) by circles and a solid line. We note that 
the positive sense of rotation 66 defined as shown in Fig.Ure- 
sulted in the negative (top-to-bottom) shift of the exit beam in 
the observation plane (see Supplementary material). The exit 


beam shifted to the bottommost position represents the up¬ 
stream portion of the first crystal. According to the diagrams 
shown in Fig. |2jc) for a convex shape of the crystal surface 
the direction of the beam shift should be opposite. Thus, the 
actual shape of the crystal in the experiment was concave. The 
resulting angular profile is plotted in FigureEIf)- It represents 
the derivative of the crystal surface profile (i.e., the slope er¬ 
ror). The surface profile in our analysis is considered with 
respect to the inward normal to the crystal surface. 

The result of FEA correctly predicts the shape yet underes¬ 
timates the measured profile. In an attempt to analyze the ori¬ 
gin of this discrepancy we established that the temperature of 
the base of the first crystal during the experiment was r; 82 K 
which agreed with the result of FEA. The maximum increase 
in the temperature with respect to this base level predicted 
by FEA was 10 K. It was concluded that the discrepancy be¬ 
tween the angular profiles can not be attributed to the result¬ 
ing change in the d-spacing (< 1/rrad) induced by such small 
temperature deviation at this low temperature. The discrep¬ 
ancy could be attributed to a residual strain due to mounting of 
the crystal or thermal destabilization of strain-free mounting 
conditions due to differential thermal expansion upon cooling 
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FIG. 5: Rocking curve topographs of the exit beam of a cryo-cooled double-crystal monochromator at an undulator beamline of the Advanced 
Photon Source (20-ID) for the case of moderate heat load (i5o^^ = 11 keV, undulator deflection parameter K = 0.75, total absorbed power P = 
93.5 W); a map of the normalized peak intensity of the local rocking curves (a), a map of the curve width (b) and a map of the peak position (c). 
Each map is sliced as shown by the vertical dashed line. The resulting corresponding distributions along these lines are those of peak intensity 
(d), curve width (e) and peak position (f). The theoretical value of the rocking curve width for perfect crystals in the double-crystal geometry 
(2.2 /irad) is shown by the dashed line in (e). The dashed line in (f) represents the reference zero level. A portion of the observation plane is 
shown which includes the central radiation cone. The total observation area is limited by an entrance aperture of 2.4 x 1.2 mm^ upstream of 
the monochromator. 


to the cryogenic temperatures. Figure|6tb) shows the angular 
prohle measured at low heat load (open gap of the undula¬ 
tor corresponding to a photon energy of 13 keV) compared to 
FEA slope error prohle calculated for this condition (absorbed 
power 24 W). The result of FEA suggests that the heat-load- 
induced slope error can be neglected. Thus, the measured 
angular prohle represents a reference level due to residual 
mounting-induced slope error which can be subtracted from 
the result of measurement at the moderate heat load (Eig- 
ure HJa)) to yield the effective heat-load-induced slope error. 
This difference shown in Eigure|6|c) by circles and the solid 
line is compared with the difference between the correspond¬ 
ing EEA-calculated prohles (dashed green line). As seen from 
the hgure such renormalization procedure yields a good agree¬ 
ment except at the downstream most portion of the beam foot¬ 
print where an asymmetric deviation from the EEA slope error 
is observed. This asymmetry could be attributed to imperfect 
alignment of the entrance aperture. Such a likely conclusion 
is supported by the fact that only the bottom beam-clipping 
edge of the aperture was observed in the experiment. 


VI. CONCLUSIONS 

In conclusion, x-ray diffraction imaging method was ap¬ 
plied for rapid in-situ diagnostics of double-crystal monochro¬ 
mators. The method is based on sequential x-ray diffraction 
topography of high-order Bragg reflections present in the exit 
beam of double-crystal monochromators at synchrotron ra¬ 
diation facilities. The method provides visualization of the 
working crystal region and its properties under the heat load 
of intense synchrotron radiation. The state of the thermally 
distorted first crystal is monitored using local rocking curves 
with a spatial resolution which can be as low as r; 10 fim. Pa¬ 
rameters of local rocking curves are mapped across the obser¬ 
vation region and can be interpreted as a map of the beam foot¬ 
print (reflected intensity), local crystal quality (curve width) 
and wavefront distortion (peak position). In the case of cryo- 
cooled working crystals the latter can be attributed to local 
misorientation of the reflecting crystal planes which can be 
compared to heat-load-induced slope errors derived from h- 
nite element analysis. The characterization procedure at any 
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FIG. 6: Angular profiles along the center line of the beam footprint on the first crystal. The results of finite element analysis are shown by 
green dashed line and the measured angular profiles are shown by circles and solid lines. The angular profiles represent the concave shape of 
the crystal. The peak in the measured profiles is due to clipping of the reflected intensity by the bottom edge of the entrance aperture. The peak 
was used for alignment purposes, (a) the first undulator harmonic is tuned to a photon energy of 11 keV; (b) the first undulator harmonic is 
tuned to a photon energy of 13 keV. (c) the experimental angular profile of (b) is taken as a reference level and subtracted from the experimental 
profile of the case (a). The resulting difference curve (red circles, red solid line) is compared with the difference between the corresponding 
results of FEA (green dashed line). 


given heat load condition involves a single angular scan of the 
second monochromator crystal and data analysis using a soft¬ 
ware code. The procedure can be easily automated. Rapid 
complete evaluation of the double-crystal monochromators 
provided by the method can facilitate development of adap¬ 
tive systems where the state of the crystal is manipulated to 
minimize heat-load-induced wavefront distortion of the exit 
beam. 
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